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Abstract 

The masses of octet baryons are calculated by the method of QCD sum 
rules. Using generalized interpolating fields, three independent sets of QCD 
sum rules are derived which allow the extraction of low-lying N* states with 
spin-parity 1/2+, 1/2- and 3/2- in both the non-strange and strange channels. 
The predictive ability of the sum rules is examined by a Monte-Carlo based 
analysis procedure in which the three phenomenological parameters (mass, 
coupling, threshold) are treated as free parameters simultaneously. Realistic 
uncertainties in these parameters are obtained by simultaneously exploring 
all uncertainties in the QCD input parameters. Those sum rules with good 
predictive power are identified and their predictions are compared with ex- 
periment where available. 
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I. INTRODUCTION 



The QCD sum rules method [jTl is a time-honored approach in reveahng a direct con- 
nection between hadron phenomenology and QCD vacuum structure via a few universal 
parameters called vacuum condensates (vacuum expectation values of QCD local opera- 
tors). The most important are the quark condensate, the mixed condensate and the gluon 
condensate. The method has been successfully applied to a variety of problems to gain a 
non-perturbative understanding into the structure of hadrons (for a review of the early work, 
see Ref. [|]), and continues to be an active field 0. 

The use of QCD sum rules to study baryon masses started not long after the method was 
introduced. Early works include Ref. i^i which a variety of baryons were studied 

in the approach. The success of this somewhat unconventional method in predicting hadron 
properties has historically given us a unique perspective on hadron structure, compared to 
other conventional views, like the phenomenologically successful constituent quark model. 

In recent years, there is increasing desire to understand the structure of baryons directly 
from QCD, the fundamental theory of the strong interaction. In particular, increasing focus 
is put on N* physics, the study of excited states baryons. After all, most of the observed 
baryon spectrum under about 2 GeV can be considered as some excitation of the ground- 
state nucleon. On the experimental side, the field is fueled by the high-quality data coming 
from JLab and other accelerators. QCD is very difficult to solve in the low-energy regime. 
While lattice QCD, which solves QCD with controlled systematic errors on a discrete space- 
time lattice by numerical simulations, holds the key to ultimately unravel hadron structure, 
its results are still hard to obtain. One the other hand, the QCD sum rule method can serve 
as a complement to lattice QCD based on the qualifications that it has close contact with 
the basic premise of QCD, has minimal model dependence with well- understood limitations 
inherent in the operator product expansion (OPE). One advantage of the method is that 
it provides a physically transparent view of how hadrons arise from QCD via the vacuum 
condensates. One can trace back, term by term, what operators and how these operators 
conspire to give rise to the mass of a baryon, and mass splittings among baryons. An example 
is the so-called loffe formula @], ~ (— Svr^ < qq >)^/^ = 0.97 GeV, which is an over- 
simplified result by keeping only the leading order contributions in the nucleon QCD sum 
rules. It nonetheless gives a nucleon mass not too far from the observed value and a direct 
link to the quark condensate, the order parameter of spontaneous chiral symmetry breaking 
of QCD. In general, the reliability of such predictions should be checked by including all 
terms and by careful analysis. But the insight is valuable in understanding hadron structure. 

To make the QCD sum rules method a useful tool in probing N* structure, a key question 
is how much predictive power the method has. Within the limitations of the method, the 
question is sum-rule dependent because all QCD sum rules are not equal, as the examples will 
show in this work. It depends on the type of correlation functions and the tensor structure 
from which a sum rule is constructed because the OPE usually has different convergence 
properties at each structure. The main goal of this work is to examine the predictive ability 
of the QCD sum rule approach in its standard implementation, using the octet baryons as 
examples. A similar study has been done for the decuplet baryons In our view, such a 
study is useful in a number of ways. 

First, we find the existing works on these baryons are incomplete and sometimes incon- 
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sistent in a number of areas, such as the interpolating fields, number of terms included in 
the OPE, the anomalous dimension corrections, the continuum model, and the treatment 
of strange baryons. A perusal of these works sometimes reveals discrepancies among the 
QCD sum rules where comparisons are possible. Therefore, before we move to the issue 
of predictive power, we need a set of sum rules as complete as possible. To this end, we 
first derive a complete set of QCD sum rules in order to independently verify and expand 
upon existing results. We use generalized interpolating fields and three types of correlation 
functions among them. We study not only the nucleon, but also the strange members of the 
octet family. Some of the sum rules in this work are new. The result is a complete, updated 
set of QCD sum rules under one roof for octet baryons in the standard implementation of 
the approach (11 sum rules for each of the iV, S, H and A channels) which encompass the 
existing sum rules as special cases. We consistently include operators up to dimension 8, first 
order strange quark mass corrections, fiavor symmetry breaking of the strange quark con- 
densates, anomalous dimension corrections, possible factorization violation of the four-quark 
condensate, and a pole plus continuum model. 

Second, it concerns how to extract a physical quantity from sum rules and how reliable it 
is. Usually there is considerable freedom in how the analysis is done regarding the Borel mass 
window, the modeling of the continuum, and the treatment of the fit parameters, each of 
which has some impact on the outcome. In the past, often an error is assigned to a quantity 
without the support of rigorous error analysis. A quantitative tool is needed to treat these 
degrees of freedom as complete as possible. This was first addressed by Leinweber [1^ who 
introduced a Monte-Carlo based procedure that incorporates all uncertainties in the QCD 
input parameters simultaneously, and maps them into uncertainties in the phenomenological 
parameters, with careful regard to OPE convergence and ground state dominance. We will 
apply this procedure in the current work. 

Third, the study will aid the application of the method to matrix elements of N* (three- 
point functions), such as magnetic moments, transition moments, axial charges, tensor 
charges, etc., because they utilize the two-point function as normalization. 

The paper is organized as follows. In section ||, the interpolating fields are introduced, 
followed by a complete listing of all the derived sum rules divided into three groups. Sec- 
tion m discusses the analysis procedure and all the results. Section gives the concluding 
remarks. 



II. METHOD 



The basic idea behind the QCD sum rules approach is the ansatz of duality. That is, it 
is possible to simultaneously describe a hadron as quarks propagating in the QCD vacuum, 
and as a phenomenological field with the appropriate quantum numbers. In this way, a 
connection can be established between a description in terms of hadronic degrees of freedom 
and one based on the underlying quark and gluon degrees of freedom governed by QCD. 

Hadron masses can be extracted from the time-ordered two-point correlation function in 
the QCD vacuum 

U{p) = i j d^x e*^'-^(0 I T{ r]{x) 7^(0) } | 0), (1) 
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where rj is the interpolating field (or hadron current) with the quantum numbers of the 
baryon under consideration. 

Assuming SU(2) isospin symmetry in the u and d quarks, we consider the most general 
current of spin 1/2 and isospin 1/2 for the nucleon: 



Ni 



where 



and 



V^ix) 



abc 



(u''^{x)Cd\x))-i^u\x). 



(2) 



(3) 



(4) 



Here C is the charge conjugation operator, the superscript T means transpose, and eabc makes 
it color- singlet. The real parameter (3 allows for the mixture of the two independent currents. 
The choice advocated by loffe @] and often used in QCD sum rules studies corresponds to 
j3 = —1. In this work, the freedom to vary j3 is exploited to achieve maximal overlap with 
the state in question for a particular sum rule. It proves an effective tool to saturate a sum 
rule either by positive- or by negative-parity states. 

In the following we write down the interpolating fields for the other members of the 
octet, omitting the explicit x dependence. For S we consider 



and for H 



-2e 



abc 



-2e 



abc 



aT/ 



(5) 



(6) 



For the A, there is the possibility of octet or fiavor-singlet quantum numbers. We consider 
three tj^es of interpolating fields for the A: the octet A, 



2y - e"^^{ [2{u^^C-f5d^)s' + {u''^C-i^s^)d^ - {d''^C-f5S^)u' 

+P [2(m'^^Cc/'')75s'= + (m"^Cs'')75c/" - (c/"^Cs^)75m"] } , 

the fiavor-singlet A, 

rj^' = 2^e"^"{ [(u'^'^C75rf^)s" - (M"^C75S^)ci" + (rf"^C75s'')u^' 

+(3 [{u^'^Cd'h^s'' - {u'^^Cs'h^d' + {d'^^Cs'h^u'] } , 
and what we refer to as the common A, 



(7) 



(8) 



2^ ^ e"'" { [(m''^C75S*)c/" - (rf'^^C75S*)M'=] + (3 [{u^^Cs'')-f5d'' - {d^^Cs^h^u"] } . (9) 
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which consists of terms common to r]^" and r]^" . Since SU(3) flavor is broken by the strange 
quark, it may be interesting to investigate such an interpolating field which does not impose 
a flavor symmetry of the quarks composing A. 

The normalization in the above octet currents (excluding rj^" and rj^") is chosen so that 
the correlation functions at the quark level coincide with each other under SU(3) flavor 
symmetry. This provides a simple check of the calculation. 

For spin-3/2 non-strange N* with isospin 1/2, we consider: 



^abc 



(10) 



It has been known that the ground-state nucleon couples to this interpolating field via its 
spin 1/2 component. As for the strange members, we consider for the E: 



„S abc 



^abc 



for the 5: 

^3/2,^ = e 
and again for the three types of A: 



^abc 



(11) 



(12) 



(13) 



and 



'l3/2,ix 



(14) 



(15) 



This completes the fist of all the interpolating fields considered in this work. 

A baryon interpolating field excites both the ground state and the excited states with 
both parities. The ability of an interpolating field to annihilate a baryon into the QCD 
vacuum is described by a phenomenological parameter (often referred to as current coupling 
or pole residue). It is defined for the spin-1/2 interpolating field by 



(0 1 771/2 1 1/2+ps) = A1/2+ w(p, s), 
for a positive-parity state with momentum p and spin s, and 

(0 1 771/2 I l/2~ps) = Ai/2- 75 u{v, s). 



(16) 



(17) 
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for a negative-parity state. Here u is the Dirac spinor with the normahzation u{p, s) u{p, s) = 
2Mb, where Mb is the baryon mass. 

For the spin-3/2 interpolating field, it is defined by 

(0 I V3/2,t, I 3/2+ps) = A3/2+ u^{p, s), (18) 

for a positive-parity state, and 

(0 I r73/2,M I S/2~ps) = A3/2- 75 u^,{p, s), (19) 

for a negative-parity state. Here is a spin-vector in the Rarita-Schwinger representation 
with the normalization Ufj,{p, s) Uf^^p, s) = —2Mb. In addition, the coupling of a spin-3/2 
current with a spin-1/2 state is defined by 

(0|773/2,^|l/2» = «i/2+ (^ + 7m) 75w(p,s), (20) 

for a positive-parity state, and 

(0 I r73/2,M I l/2-ps) = ai/2- + 7,.) u{p, s), (21) 

for a negative-parity state. Note that A, a are two different parameters coupling to spin-1/2 
states. Knowledge of these parameters is helpful for three-point calculations in the QCD 
sum rule approach, since it enters as normalization. 

The QCD sum rules are derived using standard implementation of the approach. On 
the QCD side (conventionally referred to as the LHS), the correlator in (|l]) is evaluated 
using the Operator-Product-Expansion (OPE). On the phenomenological (or the RHS), 
the correlator is populated by intermediate states. The Borel transform is applied to the 
momentum-space correlator to simultaneously provide exponential suppression of excited 
states and factorial suppression of higher-dimension operators. Furthermore a pole plus 
continuum model is usually adopted to reduce the number of phenomenological parameters 
to three: the ground-state pole mass [Mb), the current coupling (A^) and the continuum 
threshold {w). By numerically matching both sides of the sum rule over a reasonable window 
in the Borel mass parameter, these parameters are extracted as predictions of the method. 
Whether a sum rule is saturated by a positive-parity state or a negative-parity state can be 
determined by examining the leading-order terms in the LHS and the sign of RHS. 

The standard treatment of the QCD side proceeds by inserting the explicit forms of the 
interpolating fields into the two-point function of Eq. (|l|) and contracting out pairs of time- 
ordered quark-field operators which are the fully-interacting quark propagators of QCD. 
The result is a so-called master formula expressed in terms of the quark propagators. For 
completeness, all of the master formulas used in this work are given in the Appendix. 

The fully-interacting quark propagator can be derived from the Wick's theorem for the 
time-ordered product of two fermion fields under the usual assumption of vacuum saturation 
of the intermediate states of composite operators. The result, which includes operators up 
to dimension 8 and terms linear in the quark mass, is explicitly given by 
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283" " 
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A" 



(22) 



Since we deal with two independent interpolating fields for spin-1/2 and spin-3/2, we 
consider three possible types of correlations: the correlator of generalized spin-1/2 currents 
rjii2 and f/1/2, the mixed correlator of generahzed spin-1/2 current 771/2,/^ = 7^75^1/2 and 
spin-3/2 current f]3/2,u, and the correlator of spin-3/2 currents ^73/2,/^ and fjz/2,v From the 3 
types of correlators, 11 independent sum rules can be constructed which contain information 
on the lowing-lying states in the six channels (A^, S, S, Aq, A^, Ac) and in the spin-l/2± and 
spin-3/2± sectors. The resulting sum rules are given in the following sections in sufficient 
detail so that they can be checked by other researchers. 

The quark propagator is then substituted into the master formula to give the LHS of the 
sum rule. This tedious work can be organized by the 22 diagrams in Fig. ^ and Fig. ^ and 
is aided by the computer-algebra package REDUCE. 



A. QCD Sum Rules: Spin-1/2 Correlator 

In this section we give the QCD sum rules as derived from the correlator of general 
spin-1/2 interpolating fields 771/2 and f/i/2- 

The phenomenological side of the correlation function can be written as 

n(p) = ^^^(p±MB) + ---, (23) 

where the dots denote excited state contributions. The caret notation denotes p = p"7q,. 
The plus/minus sign corresponds to positive/negative parity. 

The correlation function has two tensor structures, p and 1, from which two independent 
QCD sum rules can be constructed for each member of the octet family. It turns out that 
the QCD sum rules from the structure p contain dimension-even condensates only, so we 
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will refer such QCD sum rules as chiral-even. Similarly, we refer to those from the structure 
1 as chiral-odd, because they involve dimension-odd condensates only. This classification 
scheme will be used throughout this work. 

The chiral-even sum rules coupling to spin- 1/2 states are (the tensor structure is indicated 
on the left) 

p : ci L-^/^ E2 M^ + C2b L-^'^ Eq + C3 a L'^/^ Eo 
+C4 mla L-^^'^' + C5 K,a^ L''^ + Cg my L-^'^' 
= A?/2+e 1/2+/ +A?/2_e v^-^ , (24) 
where the coefficients are given by 

N: ci = ^(5-h2/3 + 5/32), 
02 = ^(5 + 2/3 + 5/3^), 

C3 = 0, 

C4 = 0, 

C6 = -^[13-2/3-ll/32], (25) 

E: ci = ^(5 + 2/3 + 5/32), 
02 = ^(5 + 2/3 + 5/3^), 
C3 = ^[(12 - 5/,) - 2/,/3 - (12 + 5f,)(3% 

C4 = -^[(4/. + 21 + 18t(M2)) + 4/,/3 + (4/, - 21 - 18t{M^))i3% 
C5 = ^[(6/. + l)-2/3-(6/,-l)/32], 

C6 = -^[(12/3 + 1) - 2/3 - (I2/3 - l)/3^], (26) 

S: 01 = ^(5 + 2/3 + 5/3^), 
C2 = ^(5 + 2/3 + 5/32), 
C3 = ^[(2-/,)-2/,/3-(2 + /,)/32], 

C4 = -^[(15 - /. + mM^)) - 10/,/3 - (15 + /, + 18t{M^))(3% 
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C6 = -hifs + 12) - 2/,/3 + (/, - 12) P% (27) 



K- ci = — (5 + 2/3 + 5/3'), 
02 = ^(5 + 2/3 + 5/3^), 

C3 = ^[(20 - 15/,) - (16 + 6/,)/3 - (4 + 15f,)f3% 

C4 = ^[(4/. - 5 - 6t{M')) + (4 + 4/,)/3 + (4/, + 1 + 6t{M'))P% 

C5 = ^[(10/, + 11) + (2 - 8/,)/3 - (2/, + 13)P% 



C6 = ^[(-16/. - 23) + (8/3 - 2)/3 + (8/, + 25)(3% (28) 

A,: ci = ^(l-2/3 + /32), 
C2 = ^(l-2/3 + /32), 

C3 = ^[(-4 - 3/,) + (8 + 6/,)/3 - (4 + 3/,)/3^], 
C4 = 3^(1- 2/3 + /32), 

C5 = ^[(-1 - 2/,) + (2 + 4/,)/3 - (1 + 2/,)/32], 

C6 = ^[(1 + 2/.) - (4/, + 2)/3 + (1 + 2fs)f3% (29) 

Ac: ci = -^(3-2/3 + 3/32), 
Id 

C2 = ^(3-2/3 + 3/32), 

C3-^[(4-3/,) + 2/,/3-(4 + 3/3)/3^], 

C4 = ^[(-3 - 6i(M2)) + 4/,/3 + (3 + 6i(M2))/32], 

C5 = i[(2/,-l) + 2/3-(2/, + l)/32], 

C6 = ^[(-4/. + 1) - 2/3 + (4/, + 1)P% (30) 

Herein after, the condensates and couplings are re-scaled as a = — (27r)^ {uu), b = {g1 G^), 
{uQcCy ■ Gu) — —ml{uu), A = (27r)^A, a — (27r)^Q;. The ratio fs = {ss)/{uu) — 
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{sQcCr ■ Gs) I {uQcCr ■ Gu) accounts for the flavor symmetry breaking of the strange quark. We 
retain only terms linear in the strange quark mass and set m„ = = 0. The four-quark 
condensate is not well-known and we use the factorization approximation (uuuu) = ^^{uu)'^ 
and investigate its possible violation via the parameter Ky. The anomalous dimension correc- 
tions of the various operators are taken into account via the factors LP' — [as{iJ^) / as{M'^)\' — 
|^ln(M^/AQ(^^)/ln(/i^/AQ(^^) , where 7 is the appropriate anomalous dimension, ji — 500 
McV is the renormalization scale, and Aqcd is the QCD scale parameter whose value will 
be given below. The function t{M^) = In ^ — 'Jem with Jem ~ 0.577 the Euler-Mascheroni 
constant. As usual, the excited state contributions are modeled using terms on the OPE 
side surviving 00 under the assumption of duality, and are represented by the factors 

En{x) = 1 — e^^X^^x"/??-! on the LHS with x — w%/M^ and wb an effective continuum 
threshold. Note that wb is in general different from sum rule to sum rule and we will treat 
it as a free parameter in the Monte-Carlo analysis. 

A trivial check of the calculation is provided by the fact that the QCD sum rules for E 
should reduce to those for the nucleon by replacing the s-quark with the d-quark, which is 
equivalent to setting nis — and fs — 1- This check is performed throughout the work. 

The chiral-odd sum rules couphng to spin- 1/2 states at the structure 1 are 

1 : ci L"*/^ E2 M^ + C2a Ei + C3 m^a L"^^/^^ Eq 
+C4 TUg b L~^f^ Eq -\- ah + cq nis Kt,a^ 

= X%2, My,. e-"v../-^ - A?/,_ My,- e-"v.-/-^ (31) 
where the coefficients are given by 

N: ci = 0, C2 = ^[7 - 2/5 - C3 = - c, = 0, 

C5 = ^[19 + 10/5 - 29/3'], C6 = 0, (32) 

S : c, = \{l- C2 = ^[(6 + /.) - 2/,/3 - (6 - fs)P% c, ^ - P'), 

C6 = ^[(5 - 3/.) + 2/? + (5 + 3/,)/?2], (33) 

5 : = ^(1 - /32), C2 = ^[(6/, + 1) - 2/3 - (6/, - 1)/?^)], c, = -^(1 - p'), 
C4 = ^(1 - /3'), C5 = ^[(24/. - 5) + 10/3 - (24/, + b)(3\ 
C6 = ^[(3-/.) + 2/3-h(3 + /,)/32], (34) 



A„ : ci = -(11 + 2/3 - 13/3^), = -[(10 + 11/,) + (-8 + 2/,)/3 - (2 + IZfsWi 
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C3 = - 2/.) + (1 + 2f,)P'), C4 = ^(13 - 2/3 - 11/52), 

C5 = 8^[(4 + 53/.) + (40 - 10f,)P - (44 + 43/,)/?^], 

C6 = ^[(15 - 5/,) + (6 + 4/,)/5 + (15 + fs)P% (35) 



A. : Cl = -(-1 + 2/3 - /32), = -[(-2 - /,) + (4 + 2/,)/3 - (2 + /,)/32)], 
C3 = 0c4 = ^(l-2/3 + /32), 

C5 = |^[(10 + 5/,) - (20 + 10fs)P + (10 + 5fs)P% 

C6 = ^[(3 + fs) - (6 + 2/.)/? + (3 + /,)/?2], (36) 



Ac: 



Cl = 


^(-1 + 2/5- 




C3 = 




C4 = ^(l-2/3 + /32), 


C5 = 




)-10/,/?+(-8 + 5/,)/32], 


C6 = 


^[(3-/s)- 


2/3 + (3 + /,)/32], 



(37) 

Note that on the phenomenological side of the sum rules here and below, we explicitly retain 
the two lowest states, one having positive-parity and one having negative-parity, because it 
is not known apriori whether a sum rule is saturated by a positive or negative state. Only 
one term is used in the analysis once it is determined that it is the dominant one. 

B. QCD Sum Rules: Mixed Correlator 

In this section we give the QCD sum rules as derived from the mixed correlator of 
generalized spin-1/2 current ^71/2,/^ = 7^75^1/2 and spin-3/2 current f]3/2,u- 
The phenomenological side of the correlation function can be written as 

n^v(p) = p2^^^2 {tmTi^P ± Mb1^,1v T -^^i^P'^P + 27mP^| + • • • , (38) 

where the dots denote excited state contributions. The upper/lower sign corresponds to 
positive/negative parity. There are four tensor structures from which four QCD sum rules 
can be constructed for each member. It turns out that the sum rules from the structures 
lixPvP and 7^p,^ are identical, leading to three independent sum rule sets. 
The chiral-even sum rules coupling to spin-1/2 states are 
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7^7^P : cirrisa L-^/^^ Eo + cs L^^/^^ + C3 m^a L-^^/^^ + C4 wi^a^ L^/^ ^ 
= -A1/2+ ay 2+ e 1/2+ ' - Ai/2- ay 2- e v^-' , (39) 
where the coefficients are given by 

N: ci = 0, C2 = -2, C3 = 0, C4 = -^[27 - 3/?], (40) 

S : ci = C2 = -^[3(1 + /.) + (1 - 

C3 = ^[(21 + 2/,) - 2(1 + /3) i(M2) - (7 + 2/,)/?], 

C4 = ^[(13 + 14/,) - 3(2/, - 1)/?], (41) 

S : ci = -^[3(2/, + 1) + (2/, - l)/3], C2 = -^[3(1 + /,) + (/, - l)/3], 
C3 = ^[(7 + 24/,) - 2(1 + /3)(1 + 4/,) t{M^) - (5 - 8/,)/?], 
C4 = ^[(14 + 13/,)-3(2-/,)/3], (42) 

Ao : ci = C2 = -^[(1 + 5/,) - (1 - /,)/3], 

C3 = ^[(17-6t(M^)) + (5-6t(M^))/3], 

^ l^t^^ + ^^^^^ + ^"^ + ^^^^ 

A. : ci = C2 = -^[2(1 - /,) + 2(/, - l)/3], 

C3 = ^(-l + 5/?), 

C4 = ^[5(1 -/,) + 5(-l + /,)/?], (44) 
54 

Ac : ci = C2 = -^[(1 + /,) + (-1 + /,)/3], 

C3 = ^[(5-2t(M^)) + (5-2t(M2))/5], 

C4 = -^[(5 + 4/,) + (-5 + 4/,)/3], (45) 

An interesting feature of this chiral-even sum ruie is that there is oniy a reiatively smail term 
of the order rngM'^ on the OPE to model the continuum for the strange baryons, while no 
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such term for the nucleon. One consequence is that the continuum contribution is zero for 
the nucleon and small for the strange baryons, suggesting a weak coupling between excited 
states and OPE. It makes it difficult to select a valid Borel region independently to fit the 
sum rule. 

The chiral-odd sum rules couphng to spin- 1/2 states are 
7,.7z. : ci TUs M^ + C2a L^^^'^ Ei + C3 ab L^/^? ^ rUsKy L^/^^ 

= -A1/2+ ai/2+ M1/2+ e + Ai/2- ai/2- Mi/2_ e 1/2-' , (46) 

where the coefficients are given by 

TV : ci = 0, C2 = -1, C3 = -^[3 + 9/3], C4 = 0, (47) 
S : ci = C2 = -^[3(1 + /,) - (1 - /,)/3], 

4 D 

C3 = + (4 + 5/.)/3], C4 = -^(3 - /?), (48) 

S : ci = C2 = -^[3(1 + /,) + (1 - /,)/3], 

4 D 

C3 = -^^[3 + (5 + 4/,)/?], C4 = -^[3(2 + /,) + (2 - /,)/?], (49) 

Ao : ci = -i^, C2 = -^[(5 + /,) + (1 - A)/?], 

C3 = -^[(4 - /.) + (8 + /,)/?], C4 = -^(5 + /?), (50) 

A. : ci = -i^, C2 = + /.) + (1 - /.)/3], 

C3 = -2^[(1 - /.) + (-1 + /.)/3], C4 = -|(-1 + i5\ (51) 

Ac : ci = C2 = -^[(1 + /.) + (1 - /.)/?], 

C3 = -Y^[(2 - /.) + (2 + /,)/3], C4 = -^(1 + /3), (52) 

The chiral-odd sum rules coupling to spin-1/2 states are 

l^VvV ■ ci m, + C2 a L^/^t ^2 ^ ^2^ ^-2/9 

+C4 m, 6 + C5 a 6 L^/^^ ^ + ce m,«:,a^ L^/^^ ^ 

_ A1/2+ ai/2+ -m2 /m2 _ Ai/2- ai/2- -m^ /m^ 



- "^"^ e"^ 1/2-'^" , (53) 



M1/2+ Ml/2 
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where the coefficients are given by 

N: ci = 0, C2 = l, C3 = -^(3-/3), C4 = 0, C5 = --^[1 + 3/3], C6 = 0, (54) 

o 4o 

S : ci = C2 = ^[3(1 + /,) - (1 - /,)/?], C3 = -^(5 + 4/. - 3/5), 

C4 = -^(1 + 3/3), c5 = -^[3/, + (4 + 5/,)/3], c6 = -^(3-/5), (55) 

S : ci = C2 = ^[3(1 + /,) + (1 - /,)/?], C3 = -^(4 + 5/, - 3/,/?), 

^4 = 7^(1-3/5), C5 = --^[3 + (5 + 4/,)/3], C6 = -^[3(2 + /,) + (2 -/,)/?], (56) 
9o 144 o 



A. 



ci = C2 = ^[(5 + fs) + (1 - /.)/?], C3 = -^[(7 + 2/,) - (1 + 2/, 

= 2^^^ - /5)' C5 = -^[(4 - /.) + (8 + /,)/?], C6 = -^(5 + f3), 



(57) 



A. 



Cl 



12 ' 



C2 = -[(-l + /s) + (l -/.)/?], C3 



"4 = Iii(^-^)' "5 



216 



[(l-/,) + (-l + /. 



-[(-l + /,) + (l -/,)/?], 

5 

C6 = -f (-1+/5), (58) 



Ac ■■ Cl 



1-/5 



C2 = -[(l + /s) + (l-/s)/3], C3 



1 

"24 



[(1 + 2/,) + (1-2/, 



^[(2 - /,) + (2 + /,)/5], C6 = -|(1 + /5), 



(59) 



Comparing with existing calculations for the mixed correlator, the coefficients in the nucleon 



channel in the above three sum rules agree with those in Ref. The coefficients in the 
strange baryon channels are new. This turns out to be a very favorable sum rule among the 
11 sum rules in predicting the masses of spin-1/2 and positive-parity baryons. 



C. QCD Sum Rules: Spin-3/2 Correlator 

In this section we give the QCD sum rules as derived from the correlator of spin-3/2 
currents 773/2,^ and fj3/2,u- 

The phenomenological side of the correlation function can be written as 



n^.(p) = A 



3/2 



1 / N 2 PaPuP 1 2 p^Pi, 

-g^iuP + i^ilf^Pu + luP^,) + -^-J^ T Ms/2g^,^ ± -M3/27m7/^ ± ^- 



3/2 



3M 



3/2 



+ ai/2 
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4(7/.^!. + luPp,) + j^Pt^PuP + Ml/2 It^lu + j;j-^Pt^P- 



1/2 



Ml 



1/2 



+ 



(60) 
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where the dots denote excited state contributions. The upper/lower sign corresponds to 
positive/negative parity. The tensor structures are chosen because they are orthogonal to 
each other. Six independent sum rules can be constructed from the correlator for each 
member: two for spin-3/2 states only, two for spin- 1/2 states only, and two for both. 
The chiral-even sum rules coupling to spin-3/2 states are 

+ C3 a 

+C4 rus mla L'^^/^^ + C5 K,a^ L^^/^^ + my L^^/^^ -L 

= -A3/2-e '1^-' -\l/2+e 3/2+/ ^ (61) 

where the coefficients are given by 

24 5 28 

N : ci = — -, C2 = -, C3 = 0, C4 = 0, C5 = 16, C6 = — -, (62) 
00 6 

24 5 2 16 

S: ci = --, C2 = -, C3 = 8(/, + 4), c4 = --(ll/. + 28), C5 = y (4/, - 1), 

C6 = -f (4/.-1), (63) 

24 5 2 16 

S : ci = -y, C2 = -, C3 = 32, C4 = -3(5/5 + 28), cg = y/8(4 - fs), 

C6 = -^/.(4-/,), (64) 



9 



5 _8 __22 

9' " 3-^^' " 

5 _ 8 _ 20 

9' ~ s-'^' ^4 - 



16 




28 


y 




~27' 


32 




56 


¥' 


C6 = 


~27' 


16 




28 


y 


C6 = - 





(65) 
(66) 



Ac : ci = --, C2 = C3 = -Z^, C4 = -2/5, C5 = C6 = — ^, (67) 

In this sum rule, states with both positive- and negative-parity contribute with the same 
sign, which means that the excited states contamination is large. Therefore it makes it 
difficult to resolve which is the dominant contribution under the single pole-plus-continuum 
model. One could retain both terms (a two-pole plus continuum model), but it increases 
the total number of phenomenological parameters to five, which is too many for a sum rule 
to resolve. As a result, no predictions can be made from this sum rule. 
The chiral-odd sum rules couphng to spin-3/2 states are 

g^, : ci rus L'^''"' + C2 a E, L'^l-"' + C3 mla L^/^^ 

-I-C4 TTT-s C5 ao 

^16/27 

-|- Cq 777.5 l^vO- 

^16/27 

= M3/2- _ H,,, M3/2. e-^3V/-^ (68) 
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where the coefficients are given by 

_2 

N : ci = 0, C2 = 16, C3 = -8, C4 = 0, 05 = — , Cq = 0, (69) 

-^Q g 2 

S : ci = -6, C2 = y (4 - /s), C3 = —(4 - /,), C4 = 2' = -^(4 - /s), 

ce = -16, (70) 

S : ci = 24, C2 = y (4/, - 1), C3 = ^(4/, - 1), c, = -2, c, = :|(4/, - 1), 

C6 = -32/,, (71) 

16 ^ 8 ^ 1 2 ^ 16 

Ao : ci = 2, C2 = —Js, ca = --/s, C4 = --, C5 = -^/s, ce = -y, (72) 

32 ^ 16 ^ 1 4 ^ 16 

As : ci = 4, C2 = yA, C3 = — ^/s, C4 = --, C5 = -^Js, cq = — J, (73) 

16^ 8^ 1 2^ 16 

Ae : Ci = 6, C2 = y/s, C3 = --/s, C4 = --, C5 = -gJs, Cq = -, (74) 

In this sum rule, the leading terms (ci,C2) on the LHS arc positive, which leads us to 
conclude that the sum rule is dominated by negative-parity states because they come in 
with a positive sign on the RHS while the positive-parity states have the opposite sign. The 
parity cancellation in the excited spectrum may lead to a clean extraction of J''^ = | states 
from this sum rule. 

The chiral-even sum rules coupling to spin- 1/2 states are 

luP. + luPt, ■ ci L^/'^ E2 + C2 6 L^/'^ + Ca m, a L^/'^ 

+C4 rris mla L-'^/^' + c, L^^/^^ + my L'^'^' 

^-al,2+e -al/2-e v^-^ , (75) 

where the coefficients are given by 

3 5 7 

N • ci = — , C2 = - — , Ca = 0, C4 = 0, C5 = -2, ce = -, (76) 

V 3 5 (3/. + 7) 2(4/, - 1) 
ci = — ' ^2 = Ca = -[fs + 4), C4 = , C5 = , 

7( i./.s - 1) .77. 
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2 


7 




54' 


4 


7 




27' 


2 


7 




18' 



3 5 (5/. + 14) 2/,(/,-4) 
c, = -,c. = --, C3 = -4, C4 = ^^, C5 = , 

C6 = (78) 
1 5 11 

As : Ci = — , C2 = - — , C3 = --fs, C4 = -/s, C5 = Cq = — , (80) 

1 5 K 1^ ^ . 

Ci = — , C2 = - — , Ca = --/s, C4 = -/s, ^5 = -g, C6 = — , (81j 

Note that the RHS is always negative, but the leading term (M^) on the LHS side is positive. 
This is an indication of a poor sum rule. In fact, no physical results have been extracted 
from this sum rule. 

The chiral-odd sum rules coupling to spin- 1/2 states are 

+ C2 mla 
+C3 rus + ab 

^16/27 

+ C5 rUs KyQ, 

^16/27 

= -q;?/2+ ^1/2+ e 1/^+^ + q;^/2- Mi/2_ e V^-' , (82) 

where the coefficients are given by 

12 5 
^ ■■ ci = -, C2 = --, C3 = 0, C4 = — , C5 = 0, (83) 

1 2 15 

S : ci = --(/, - 4), C2 = g(A - 4), Ca = —,04 = -^(/. - 4), 

C5 = -^(l + /.), (84) 

1 2 15 

S : ci = -(4/, - 1), C2 = --(4/, - 1), Ca = --, C4 = ^(4/. - 1), 

C5 = -^/s(3 + 2/,), (85) 



1 X 2 ^ 1 5 ^ 4 

Ao : ci = — C2 = -^/., Ca = - — , C4 = — C5 = --, (86) 



A 2 „ 4 „ 1 5 „ 4 

ci = — /s, C2 = Ca = - — , C4 = ^/s, C5 = --, (87) 
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Ac : ci = -/„ C2 = --/., C3 = --, C4 = — /„ C5 = --, (88) 

The positive coefficient of tfie leading term suggests that this sum rule should be saturated 
by a negative-parity state. However, in the ol:)served spectrum the state is lower than 
the I state. Since here one does not have the freedom of varying the interpolating field to 
minimize the influence of the positive-parity term, the positive-parity term must be retained 
in a sensible fit. This worsens its predictive ability. Without knowing the order of the two 
states, it is the natural position to take to keep both terms. 

The chiral-even sum rules coupling to both spin-3/2 and spin- 1/2 states are 



_ 2 -^3/2- -M|^^_/M2 2 -^3/2+ -M|^^_^/M2 

+16 ^ e-^V/^^ + 16 ^ e-^?/-/^^ (89) 
^1/2+ -'^1/2- 



where the coefficients are given by 



N: ci = ^, C2 = ^, C3 = C4 = 0, (90) 
o o 

(91) 
(92) 
(93) 
(94) 
(95) 

The chiral-odd sum rules couphng to both spin-3/2 and spin- 1/2 states are 
P„Pu ■■ ci TUs El M^ + C2a L^^^^'^ Eq + C3 m^a L^/^^ 

+C4 m, b L-'/'' + c,ab L^'l^' + ^6 k^o" L'''^' 

2 '^3/2- _/M2 2 -M2 , /M2 

O iW3/2- -J -'^'^3/2+ 

-8 ^ e-</^-/^^ + 8 ^ e-^v.-/*^^ (96) 
^V2+ ^V2- 



E : 


Cl = 


12 

= — , C2 
5 ' 


1 

^3' 


-8/„ 


C4 = 


3/., 




Cl = 


12 


1 


-16/., 


C4 = 


20 


Ao: 


Cl 


4 

= C2 
5' 


1 

^9' 




C4 = 




A,: 


Cl 


4 

= 5' 


1 

^9'"^ = 


-3/., 


C4 = 




Ae: 


Cl 


4 

= C2 
5' 


1 


3/., 


C4 = 


3/., 
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where the coefficients are given by 



N: ci = 0, C2 = -8, C3 = 8, C4 = 0, cg = -1, Ce = 0, (97) 

8 8 11 "^2 

E : ci = 2, C2 = -(/, - 4), C3 = --(/, - 4), C4 = -,05 = -(/, - 4), c, = -/„ (98) 

6 6 b 6 6 

8 8 2 1 

S : ci = -8, C2 = -3(4/^ - 1), C3 = -(4/, - 1), C4 = --, C5 = -3(4/^ - 1), 

16 

C6 = y/.(2/. - 1), (99) 

2 8 8 1 1 

Ao : ci = --, C2 = --/s, C3 = -/,, C4 = - — , C5 = -gfs, Cq = 0, (100) 

4 16 16 12 

As : ci = --, C2 = — ^/s, C3 = y/s, C4 = --, C5 = --/s, C6 = 0, (101) 

8 8 11 

Ac : ci = -2, C2 = --fs, C3 = -/s, C4 = --, C5 = --/s, C6 = 0, (102) 

3 3 6 

Note that we exphcitly keep the four lowest states: 3/2+, 3/2— , 1/2+, 1/2— . That too 
many unknown parameters are present in the two sum rules makes it impossible to isolate 
the different states in a straightforward Monte-Carlo analysis. Their usefulness only lies in 
serving as a consistency check using the parameters extracted from the other sum rules. 



III. ANALYSIS 

To extract physics information from the sum rules, we employ a Monte-Carlo-based anal- 



ysis procedure ffist introduced in |10|. The basic steps of the procedure are as follows. Given 
the uncertainties in the QCD input parameters, randomly-selected, Gaussianly-distributed 
sets are generated, from which an uncertainty distribution in the OPE can be constructed. 
Then a minimization is applied to the sum rule by adjusting the phenomeno logical fit 
parameters. This is done for each QCD parameter set, resulting in distributions for phe- 
nomenological fit parameters, from which errors are derived. Usually, 100 such configura- 
tions are sufficient for getting stable results. We generally select 500 or more sets which 
help resolve more subtle correlations among the QCD parameters and the phenomenological 
fit parameters. One distinctive advantage of the method is that the entire phase-space of 
the input QCD parameters is explored simultaneously, leading to more realistic uncertainty 
estimates in the phenomenological parameters. 

In addition, the Borel window over which the two sides of a sum rule are matched is 
determined by the following two criteria: a) OPE convergence — the highest-dimension- 
operators contribute no more than 10% to the QCD side; b) ground-state dominance — 
excited state contributions should not exceed more than 50% of the phenomenological side. 
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The former effectively establishes a lower limit, the latter an upper limit. Those sum rules 
which do not have a valid Borel window under these criteria are considered unreliable. The 
emphasis here is on exploring the QCD parameter space via Monte Carlo. The 10%-50% 
criteria are a reasonable choice that provide a basis for quantitative analysis. Reasonable 
alternatives to the criteria are automatically explored in the Monte-Carlo analysis, as the 
condensate values and the continuum threshold would change in each sample. 

The QCD input parameters and their uncertainty assignments are given as follows. The 
quark condensate in standard notation is taken as a = 0.52 ± 0.05 GeV^, corresponding 
to a central value of {uu) = —(236)^ MeV^. For the gluon condensate, early estimates 
from charmonium place it at 6 = 0.47 ± 0.2 GeV^, a value commonly used in QCD sum 
rule analysis. But more recent investigations support much larger values [|12]-[1^. Here we 
adopt 6 = 1.2 ± 0.6 GeV^ with 50% uncertainty. The mixed condensate parameter is placed 
at rriQ = 0.72 ± 0.08 GeV^. For the four-quark condensate, there are claims of significant 
violation of the factorization hypothesis ]T^|-p!^. Here we use k^, = 2 ± 1 and 1 < < 4. 
The QCD scale parameter is restricted to Aq(7£)=0.15±0.04 GeV. We find variations of 
Aqcd have little effects on the results. The strange quark mass is taken as = 0.15 ± 0.02 
GeV. The value of fg has been determined in and is given by fs = 0.83 ± 0.05 after 
converting to our notation by 7 = /s — 1. These uncertainties are assigned conservatively and 
in accordance with the state-of-the-art in the literature. While some may argue that some 
values are better known, others may find that the errors are underestimated. In any event, 
one will learn how the uncertainties in the QCD parameters are mapped into uncertainties 
in the phenomenological fit parameters. 

Before going into the numerical analysis of the sum rules, a general discussion of the 
reliability of QCD sum rules is in order. It has been argued that chirality plays an important 
role in determining the reliability of baryon QCD sum rules JlBl. One reason is that the 
contributions of positive- and negative-parity excited states partially cancel each other in 
the chiral-odd sum rules and add up in the chiral-even sum rules. An inspection of the 11 
sum rules in this work indeed exhibits this feature. In the six chiral-odd sum rules, Eq. (|3TD 
from the spin-1/2 correlator, Eq. (|46|) and Eq. ( ^31) from the mixed correlator, and Eq. (|68|), 
Eq. (0), and Eq. ( p6D from the spin-3/2 correlator, the positive- and negative-parity states 
have the opposite sign, whereas in the chiral-even sum rules, they have the same sign. As 
a result, the continuum contamination is much smaller in the chiral-odd sum rules than 
in the chiral-even ones. In addition, the chiral-even sum rules suffer from other relatively 
large uncertainties, arising from large perturbative corrections in the OPE, the factorization 
assumption of chiral-even condensates, and the onset of non-factorizable condensates in 
relatively low dimensions §]. So it is expected that the chiral-odd sum rules are generally 
more reliable than the chiral-even sum rules before a numerical analysis. 



A. Predictive Power: Three Parameter Search 



A three-parameter search is by treating all three phenomenological parameters (contin- 
uum threshold w, coupling and baryon mass Mb) as free parameters simultaneously. 
Since the QCD input parameter space is also explored simultaneously by the Monte-Carlo 
approach, a successful three-parameter search provides the best demonstration that a QCD 
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TABLE I. Predictions from 3-parameter search of the chiral-odd sum rules Eq. (^T]), at the 
structure 1 from the spin-1/2 correlator. The Borel regions determined by OPE convergence and 
pole dominance are also given. Experimental masses are taken from |17|. 



Sum 


Rule 




Region 
(GeV) 


w 
(GeV) 


(GeV^) 


Mass 
(GeV) 


Exp. 

(GeV) 




(/5 = 


-0.8) 


0.75 to 0.91 


1.24± .16 


0.70± .26 


1.02± .11 


0.938 


Si 

2 


i(3 = 


-6) 


1.00 to 1.37 


1.96± .57 


87.2± 53.2 


1.54 ± .34 


1.620 


2 ^ 


(/? = 


-6) 


1.08 to 1.22 


2.03± .68 


109.6± 84.3 


1.55± .38 


1.69(?) 




iP = 


-6) 


1.08 to 1.12 


1.72± .46 


60.3± 44.4 


1.38± .32 


1.670 


*2 


(/? = 


-6) 


1.18 to 1.24 


1.78± .20 


41. 2± 16.3 


1.43± .13 


1.405 




iP = 


-6) 


1.10 to 1.25 


l.SOib .23 


49. 2± 22.5 


1.46± .14 


1.405 



sum rule has predictive power. 

We performed a sweeping analysis of all the 11 sum rules for each baryon member. As 
expected, we found that the chiral-odd sum rules perform better than the chiral-even sum 
rules. Only three sum rules, all chiral-odd, allow a 3-parameter search. The results are 
discussed below. 

The first is Eq. (31) from the spin-1/2 correlator. Traditionally, it is a controversial sum 
rule regarding its performance. For a long time, its partner, the chiral-even Eq. (0), was 
identified as the better sum rule with the value of (3 taken as —1. It was proven otherwise 
in the Monte-Carlo analysis of [|1^]. Moreover, for the chiral-odd sum rule, the optimal 
mixing for the nucleon was determined to be /3 = —0.8 rather than —1. This result for 
the nucleon is reproduced in this work. On the other hand, in our analysis of the S, S, 
and A states, we could not find a solution when we assume the sum rule is saturated by a 
positive-parity state while tuning /3 in a wide range. Instead we found that the 3-parameter 
search is successful when we assume the sum rule is dominated by a negative-parity state 
for f3 > —2. The results for /3 = —6 are given in Table |l[ The masses and thresholds are 
quite robust, insensitive to variations in f3 in the range of —2 to —20. The only thing that 
changes with P is the coupling as it should be /? dependent. So we conclude that this sum 
rule predicts a positive-parity state in the non-strange channel with j3 = —0.8, but the same 
sum rule predicts negative-parity states in the strange channels with a large negative value 
of (3 which means a large negative component of X2 in the interpolating field in Eq. (^. 
The uncertainties were derived from 500 sets of QCD parameters by Monte-Carlo. The 
predictions compare favorably with experiments. In the case of the QCD sum rule 
predicts a mass of about 1.60 GeV. The Particle Data Group [jl^ lists two 3-star S states 
with unknown spin-parity, 1690 MeV and 1950 MeV. Our prediction favors the lower one. 
The three Lambda types give roughly the same masses within uncertainties. 

The above is an example of how a sum resolves the parity of a state. Recall that a 
baryon interpolating field couples to both positive- and negative-parity states. That is why 
two terms on the phenomenological side are explicitly kept, one of each parity. But keeping 
two terms increases the number phenomenological parameters from 3 to 5, which is too 
many in practice. Rarely a QCD rum rule has enough information to yield a successful 
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TABLE II. Predictions from 3-parameter search of the chiral-odd sum rules, Eq. (53), at the 
tensor structure "f^PuP from the mixed correlator. 



s 


um 


Rule 


Region 


w 


-^1/201/2 


Mass 


Exp. 








(GeV) 


(GeV) 


(GeV6) 


(GeV) 


(GeV) 


2^ 


(/3 


= +1.0) 


1.06 to 1.46 


1.31+ .22 


1.13+ .53 


1.06+ .11 


0.938 


2^ 


(/3 


= +1.0) 


1.12 to 1.53 


1.48+ .23 


1.62+ 0.68 


1.16 + .12 


1.193 




(/3 


= +1.0) 


1.35 to 1.80 


1.69+ .27 


2.61+ 1.20 


1.32+ .14 


1.318 




(/3 


= +1.0) 


1.27 to 1.74 


1.53+ .24 


0.66+ 0.28 


1.23+ .12 


1.116 


'^4* 


(/? 


= +1.0) 


1.28 to 1.72 


1.53+ .24 


0.66+ 0.28 


1.23+ .12 


1.116 



5-parameter search. So the standard practice is to use a single pole plus continuum model 
with 3 parameters to maintain a reasonable level of predictive power. But which of the two 
terms dominates a QCD sum rule depends on the interpolating field. Tuning (3 changes 
the makeup of the interpolating field, therefore alters the coupling strength to the states. 
The advantage of having an adjustable interpolating field is that one could vary it to have 
maximal coupling to a state of definite parity. As a consistency check, we kept both terms 
and search the couplings XI/2+, •^1/2-; threshold w simultaneously by setting the masses 
for N1/2+ and N1/2- at the experimental values. The result is that the coupling XI/2+ comes 
out several orders of magnitude greater than Xi/2-j showing that the sum rule is indeed 
dominated by the N1/2+ state. The same check shows that for the strange baryons, the 
negative-parity states dominate. 

The next chiral-odd sum rule that permits a 3-parameter search is Eq. (^3]) at the tensor 
structure "y^PuP from the mixed correlator. The results are shown in Table ^ This is 
probably the best sum rule in predicting the = states. The Borel windows are 
relatively wide. The predicted masses are reasonably close to the experimental values. And 
the results are stable, and insensitive to the value of /3 as long as /5 > 0. The sum rules 
for Aq1_,_ and A^i^_ give nearly identical results and the masses slightly overestimates the 
experimental numbers. We could not find a solution for the singlet Agi^_. The uncertainty 
on the masses from the Monte-Carlo analysis is on the order of 100 MeV. 

The other chiral-odd sum rule from the mixed correlator, Eq. (|46|) at the jfj.ju struc- 
ture, did not allow a 3-parameter search under the established fitting policy. Our analysis 
showed that the continuum contributions in these sum rules are relatively large (over 50% 
contribution), thus no valid Borel windows exist. However, if we relax the pole-dominance 
criterion and allow continuum contributions to exceed 50%, we could obtain reasonable and 

stable 3-parameter search results for = -k^ . The results will not be shown because of 

1 • • ' 

this compromise. 

Next we examine the chiral-odd sum rule of Eq. (|68|) at the structure g^^i, from the spin- 
3/2 correlator, which couples only to spin-3/2 baryon states. It allows a 3-parameter search 
and is saturated by = ^ states. One feature here is that there is no tunable parameter 
in the spin-3/2 correlator, which leaves the Monte-Carlo analysis quite straightforward. The 
results are given in Table |IT1[ If we explore the contribution from each term of the OPE, we 
can easily see by comparing E and S that the prominent difference is on the leading term, 
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TABLE III. Predictions from 3-parameter search of the chiral-odd sum rules, Eq. (pq) at the 
structure gfj,p from the spin-3/2 correlator. 

Sum Rule 

2 

S3_ 
2 

2 

"2 

* 2 

^2 



Region 
(GeV) 


w 
(GeV) 


^^3/2 

(GeV6) 


Mass 
(GeV) 


Exp. 
(GeV) 


0.95 to 1.17 


1.65± 


.24 


27.6± 11.8 


1.44± .13 


1.520 


1.29 to 1.36 


1.91± 


.25 


46.6± 20.1 


1.69 ± .14 


1.580 


1.30 to 1.39 


2.19± 


.27 


84.8± 42.9 


1.84± .16 


1.820 


1.48 to 1.55 


2.36± 


.35 


12. 7± 7.4 


2.00± .21 


1.690 


1.30 to 1.39 


2.12± 


.28 


16.8± 8.7 


1.80± .16 


1.690 


1.22 to 1.32 


2.01± 


.25 


19.8± 9.3 


1.71± .15 


1.520 



which originates from flavor-symmetry breaking via the strange quark mass m^. For S, this 
term acts to decrease the mass of S3_, while S3_ it acts to increase its mass. 

2 2 



B. Predictive Power: Multiple Sum Rule Search 

Besides the 3-parameter search for a single sum rule, we made further efforts to explore 
the predictive ability of the QCD sum rules. One possibility that is afforded by the Monte- 
Carlo method is to perform a simultaneous search of multiple sum rules, each having an 
independent Borel window and continuum threshold. The benefit is that more information 
is introduced on the OPE side and at the same time the number of parameters that are 
shared across the multiple sum rules is reduced. This may lead to increased stability and 
smaller errors in the output. In principle, it should be done with those sum rules that come 
from the same correlator. The only successful multiple sum rule search is from the sum rules 



of Eq. (^) and Eq. (|46|) from the mixed correlator, as shown in Table It is essentially 
a 5-parameter, 2-sum rule search for the = \ states. The results improve upon those 
in Table In the case of flavor- singlet A, it predicts a mass of around 2.6 GeV. So QCD 
sum rules predicts the correct ordering of flavor- singlet Lambda with | state of 1.405 GeV, 
and 1^ state of much higher mass. This relatively high value for the = |^ flavor- singlet 
Lambda is yet to be identified by experiment. The octet A and the common A interpolating 
fields produce comparable results. 



C. Two Parameter Search Results 



Absence of a 3-parameter search, the alternative is to perform a 2-parameter search 
by fixing one of the parameters to a certain value. In our analysis we use the following 
procedure. We start by fixing the baryon mass at the experimental value and search for the 
continuum threshold and coupling simultaneously. Then we fix the continuum threshold to 
the returned value and search for the mass and the coupling. To make sure the consistency 
of above two searches, a third search for the threshold and the mass is conducted by fixing 
the coupling to the returned value. Fitting results of the three parameters are accepted if 
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TABLE IV. Results from a 5-parameter, simultaneous search of the two chiral-odd sum rules 
of Eq.( |4^) at the structure '^^x'^v (labeled as 1) and Eq.( ^) at the structure 'y^p,/P (labeled as 2) 
from the mixed correlator. 



Sum Rule 


Regionl 


wl 


Region2 


w2 




-^l/2+«i/2+ 


Mass 


Exp. 






(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV^) 


(GeV) 


(GeV) 






0.69 to 1.12 


1.44 ±.15 


0.80 to 1.52 


1.35 ± 


.16 


1.07±.25 


1.02± .07 


0.938 


2^ 


(/3 = -i) 


0.75 to 1.24 


1.67 ± .25 


0.88 to 1.68 


1.58 ± 


.25 


1.72 ± .66 


1.17± .11 


1.193 




(/? = -!) 


0.97 to 1.34 


1.89 ± .32 


1.10 to 1.81 


1.80 ± 


.30 


2.70 ± 1.23 


1.32± .14 


1.318 




(/3 = -l) 


0.74 to 1.24 


1.67 ± .25 


0.88 to 1.66 


1.58 ± 


.25 


0.57 ± .22 


1.17 ± .11 


1.116 




(/? = -!) 


1.41 to 1.98 


3.46 ± .28 


1.64 to 2.45 


3.55 ± 


.39 


1.71 ± .72 


2.65 ± .21 


? 




(/?=-!) 


0.80 to 1.06 


1.53 ± .45 


0.86 to 1.21 


1.27 ± 


.08 


0.38 ± .05 


1.09 ± .03 


1.116 



above three searches yield consistent results. In the absence of such consistency, we reduce 
to giving results from just fixing the continuum threshold and searching for the mass and 
the coupling, which is conventionally done. 

Results from such two-parameter searches should be treated with caution because they 
are not truly predicted in the sense established earlier. It is no secret that for most sum rules 
it is possible to produce the desired mass by fixing the continuum threshold to a favorable 
value and adjusting the Borel window accordingly. In such cases, it is possible that the 
continuum contribution is too large, and the sum rule is analyzed deep in the perturbative 
regime, away from the cross-over region between perturbative and non-perturbative where 
the QCD sum rule is based. 

As an example, we turn our attention to the chiral-even sum rule of Eq. (|24D at the p 
structure from the spin-1/2 correlator. This sum rule has been traditionally identified as the 
more reliable of the two sum rules from the spin-1/2 correlator. However, it was shown by 
Leinweber fl^ that a native 3-parameter search with the conventional value of = — 1 gives 
a prediction of the nucleon mass around 450 MeV, which is very different from the observed 
value. Closer examination revealed that the sum rule suffers from large uncertainties in the 
OPE and a lack of a valid Borel window, so a three-parameter search is unstable. However, 
we found that a two-parameter search could produce sensible results for this sum rule, as 
shown in Table |V|. We set the beta value to the traditional value of /3 = —1. The results 
compare favorably with experiments. Again it favors a high value of around 2.6 GeV for the 

= 1^, flavor-singlet A. No solution could be found for the common A. 

Another sum rule that can produce sensible two-parameter search results is the chiral- 
odd sum rule of Eq. (|46|) at the structure 7^7^, as shown in Table t |VI| . The results have 



broad agreement with those in Table |V|. The fits are very stable and have smaller than usual 
uncertainties in the parameters. 

The chiral-even sum rule of Eq. (|7^) at the structure •y^Pu + luPfj. and the chiral-odd 
sum rule of Eq. ( ^2|) at the structure 'j^'ju + ^Qfiu in the spin-3/2 correlator couple only to 
J = ±1/2 baryon states. Numerical analysis confirms that the chiral-even one is a poor 
sum rule with nothing produced. The chiral-odd one is dominated by both positive-parity 
and negative-parity states. The Monte Carlo fit is unable to be performed without the 
positive-parity term. We tested this by fitting the couplings a^^g+j ^^1/2-^ threshold 
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TABLE V. Two-parameter search results from the chiral-even sum rules of Eq. (^) at the 
structure p from the spin-1/2 correlator. 



s 


um Rule 




Region 
(GeV) 


w 
(GeV) 


5^2 

\/2 

(GeV^) 


Mass 
(GeV) 


Exp. 

(GeV) 


Ni 

2 


, (/3 = - 


1) 


0.77 to 1.28 


1.88± .17 


1.97± .18 


0.90± .10 


0.938 


2^ 


- (/? = - 


1) 


0.79 to 1.40 


2.14± .10 


3.45± .25 


1.15 ± .17 


1.193 




. (/3 = - 


1) 


0.74 to 1.43 


2.27± .15 


4.50± .94 


1.21± .11 


1.318 




4- (/?=- 


-1) 


0.73 to 1.35 


2.07± .16 


2.76± .16 


1.08 ± .09 


1.116 


^4- 


^ (/3 = - 


-1) 


1.19 to 1.98 


3.16 ±.33 


51.5±19.0 


2.64 ± .16 





TABLE VI. Two-parameter search results from the chiral-odd sum rules of Eq. (46) at the 
structure ^^'yp from the mixed correlator. 



Sum Rule 


Region 


w 


Al/2«l/2 


Mass 


Exp. 




(GeV) 


(GeV) 


(GeV6) 


(GeV) 


(GeV) 


iVi+ (/3 = 0) 


0.68 to 1.03 


1.31± .03 


0.81± .90 


0.94 ± .02 


0.938 


Si+ (/? = -1) 


0.756 to 1.28 


1.72 ± .07 


1.81± .30 


1.18± .04 


1.193 


(/3 = -l) 


0.976 to 1.36 


1.88 ± .06 


2.55± .42 


1.31 ± .03 


1.318 


Aoi+ (/3 = -l) 


0.75 to 1.19 


1.60 ± .08 


0.49± .09 


1.11 ± .04 


1.116 


A4+ (/? = -1) 


1.42 to 1.96 


3.30 ± .39 


1.38± .67 


2.67 ± .19 




Aei+ (/?=-!) 


0.80 to 1.03 


1.45 ± .05 


0.34± .05 


1.12 ± .02 


1.116 



w simultaneously and fixing the masses of -B1/2+ and -B1/2- at experimental values. The 
two coupling constants came out very close to each other, which means the positive- and 
negative-parity states are equally important in this sum rule. We could obtain reasonable 
results for negative-parity states by fixing the positive-parity masses at their observed values. 
But the continuum contributions were found to be too large in the fits. Therefore we decided 
not to include the results here. 

The last two sum rules of Eq. (^9D and Eq. (|96D involve a mixture of | and \ states, 
which makes them not suitable for a Monte-Carlo search in their own right due to the 
numerous unknown parameters. Their utility lies in serving as a consistency check of the 
other sum rules. For example, when we fixed the parameters of the states in the chiral- 
odd sum rule of Eq. ( p6D to the values determined from the other sum rules, and searched 
for the continuum threshold and the | parameters, we obtained results that are consistent 
with those in Table |IT1[ 



IV. SUMMARY AND CONCLUSION 



We have performed a comprehensive study of all 11 QCD sum rules using general inter- 
polating fields for spin-1/2 and spin-3/2 baryon states. They are derived from three types 
of correlation functions: 2 from the spin-1/2 correlator, 3 from the mixed correlator, and 
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6 from the spin-3/2 correlator. In terms of the chirahty of the QCD operators in the sum 
rules, the breakdown is that 1+2+3=6 are chiral-odd (which contain dimension-odd QCD 
condensates only), and 1+1+3=5 are chiral-even (which contain dimension-even QCD con- 
densates only). Our results lend support to the general statement that chiral-odd sum rules 
are better than chiral-even sum rules as far as baryon two-point functions are concerned [ITB 



The sum rules contain information on the masses of A^, E, H and A states with spin- 
parity = \ and = \ ■ A Monte-Carlo based method is used to assess whether a 
given QCD sum rule has good predictive ability which we define as allowing a 3-parameter 
search in a minimization. We found that 3 sum rules, all chiral-odd, Eq. (|3T|) from the 
spin-1/2 correlator, Eq. (|53|) from the mixed correlator, and Eq. (^) from the spin-3/2 
correlator, have good predictive ability. They give predictions for the | states in Table |, 
for 1^ states in Table 0, and for | states in Table |T|. Furthermore, the Monte-Carlo 
method affords the ability to fit multiple sum rules at the same time. A multiple search 
of the two chiral-odd sum rules from the mixed correlator yields even better results for |^ 
states as given in Table |V. 

The chiral-odd sum rules from the mixed operator produce sum rules that have better 
OPE convergence which in turn lead to robust predictions. The results bode well for the 
use of mixed correlators in future studies. We also performed conventional, two-parameter 
searches by fixing one of the parameters. We can reproduce the results in the literature, and 
obtain sensible new results. But such results have certain bias built-in from the start and 
should be taken with caution. 

The predictions compare favorably with the observed values, with an uncertainty on the 
order of 100 MeV. The couplings come as by-products which are useful in the calculation of 
matrix elements because they enter as normalization. In the case of flavor- singlet A with i"*", 
the QCD sum rule approach predicts a high value of around 2.6 GeV. In the case of S, the 
QCD sum rule approach favors the state with 1690 MeV as having over the heavier state 
of 1950 MeV, both of which have unknown spin-parity assignment in the particle data table. 
Taken together, the results show that the QCD sum rule approach has predictive power for 
the masses of the excited baryons considered in this work, but that such predictive power 
must be examined on a sum rule by sum rule basis with rigorous analysis. 

One drawback in the current approach is that states with both parities contribute in the 
same sum rule for a baryon. Although sometimes one can saturate a sum rule with a certain 
parity, by using a general interpolating field with a adjustable parameter (like /3), it is not 
always feasible. For example, no information could be extracted for the |^ states because 
no sum rules from the interpolating fields considered in this work are dominated by them. It 
is desirable to perform a parity separation, as proposed in Ref. which can be achieved 
by replacing the time-ordering operator T in the correlation function in Eq. (|l]) with xq > 0, 
and constructing sum rules in the complex po-space in the rest frame (p = 0). Work that 
applies the parity projection method to the sum rules in this work and that analyzes them 
with the Monte-Carlo procedure will be reported elsewhere JT^|. 
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APPENDIX A: MASTER FORMULAS USED IN THE DERIVATIONS 

In this appendix we collect the master formulas used in the derivation of the QCD Sum 
Rules. They are functions of the quark propagators, obtained by contracting out pairs of 
time-ordered quark- field operators in the two-point correlation function in Eq. The 
same master formulas are used in lattice QCD calculations where the fully-interacting quark 
propagators are generated numerically by Monte-Carlo, instead of the analytical ones used 
here as in Eq. (|22|) . 

1. Spin-1/2 Correlators 

Using the spin-1/2 interpolating field of Eq. (0) for the nucleon (with uud quark content) 
, we find 

(r]|r{<,(x)r/,^,(0)|fi) = -4e«V''''='{ 
+Sf^,CSf^C^,Si'' + 5r'Tr(C5f ^C75^f 75) 
+/3^,Sf^,CST'''cSt'' + /?75^r'Tr(CSf ""csf 75) 
+(3SfCSf^C^,S':''i, + psf^.TriCSf'^C^.S^'') 

+P'^,SfCSf'^CSt''^, + P'^,Sf^,MCSf''cSt'') }. (Al) 

The master formula for S (with uus quark content) can be readily obtained by replacing 
the d quark with a s quark. Similarly, the master formula for S (with ssu quark content) 
can be obtained by a two-step replacement in Eq. ( |A1|) : u replaced with s, followed by d 
with u. The master formulas from the three A (with uds quark content) interpolating fields 
have more complicated expressions, for the octet A: 

{n I T{ r/f;,(x) } I ^) = -^e^^^e-'^'^'l 

+2Sf^,CSf^C^,S':'' - Sf^,CSf^C^,St'' + 5r'Tr(C5r''^C755f 75) 
+2Sf^,CSf''c^,Sf - ST'^,CST'^C^,Sf + ^r'Tr(C5r'^C755f 75) 
+2psfCS^fc^,Sf^, + 2psfCSf^C^,S'^'^, + A^Sf ^.TiiCSf^ C^,S'^') 
+2^ST'CS^fc^,Sf'^, - psfCSf^C^.Sf^, + psf^,TriCSf^C^,S^'') 
+2psfCSf'^C^,S':''^, - pST'CSf^C^.Sf^, + /?^r'75Tr(C^r'^C755f ) 
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+2p%ST'CSf^CS':''^, - P^^.SfCSf^CS'^'^, + /3%55r'75Tr(C5r'^C5f ) 
for the singlet A: 

{n I T{ T^tM^) t;M } I n) = -l^-V^'^'i 

-Sfl.CS^'^C-i.Sf - Sfl.CS'-fc^.S'^' + 5r'Tr(C5f ^C75^f 75) 
-Sf-l,CS''fc^,Sf - ST'l,CST'^C^,S'^' + 5r'Tr(C5r'^C755f 75) 
-S^'-i,CSf^C^,Sf - Sf^,CSf^C^,Sf + 5r'Tr(C^r'^C^75^f 75) 
-f5SfCS'^fc^,Sf^, - psfCSf^C^.St''^, + psf^.TiiCSf^C^.Sf) 

-psfCST''c^,Sf^, - (3ST'CSf^C^^Sf^, + (3SZ'l,T,{CSf''c^^Sf) 

-P^sST'l^CSf^CSf -P^,Sf^,CSf^CSl'' +P^,Sr'Tr{CSf^CS^''^,) 

-P'^.SfCSf^CSf^, - p'^,SfCSf^CSf-f, + f3'^,Sf-f,Tr{CSf'^CSf) }, 
and for the common A: 

(fi I T{ r/f;,(a;) } \ fi) = -2e«''^e«'^'^'{ 

-ST'l,CSf^C^,St'' + Sf Tries f'^C^.S'/^,) 
-Sf^,CSf^Cj,S'/ + SfTr{CSf^C^,Sf^,) 
-psfCSf^C^.S^''^, + /35r'75'IV(C5r'^C75^f ) 

-h^Sf^.CSf'^CS'^' + /375^r'TV(C5r'^C'5f 75) 

-P'^.ST'CSf^CSi''^, + p'^,ST'l,Tr{CSf^CSt'') 
-P%SfCSf^CSf^, + (3%Sf^,Tr{CSf^CSf) }. 
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2. Mixed Correlators 



The mixed correlators come from the generalized spin-1/2 current ?7i/2,/i = 7/^75 ^1/2 with 
rji/2 as given in Eq. (|^), and the spin-3/2 current fjz/2,u as given in Eq. (|10|). For the nucleon, 
the master formula is given by 



-(3^,Sfa,,CST'^CSi''^,a^' + (3^,ST'^.a,,TT{Sfa'''CS:^'^C) }. (A5) 

The master formulas for S and S are obtained in the same way as the spin-1/2 case described 
in the previous section. For the octet A, we have 

+27^755rV,AC'5:^'^C755f 7.^''" + 27M755rV,AC'5f ""^75^^ 7.^''" 
+47^755r'7.^pATr(755fa^^C5f^C) 

+2(3^^SfapxCS':'^CSf^,aP^ + 2(3^^Sf a^^CS^f CS"^^' ^^a^^ 
+Af3^^Sf ^,ap^i:i{Sf aP^CS':''^ C) 

+2P^,ST'a,^CS::'^CSf^,a''^ + [5^,8^ a.^CSf^ CS'^' ^^a^^ 
+f5^,ST'lucypx'TT{SfaP^CS::'^C) 

-2(3^,Sfa^,CSf^CS':''^,aP' + P^.Sf a^,CSf^CS^''^,a''' 
+f3i,Sf^,a^,TT{SfaP'CSf^C) }. (A6) 

For the singlet A, we have 

2 



-7^75 5; 

+7a<75^; 
-7^75^7' 
+7^75^7' 



,a,,Tr(755f a^^CSf C) 



30 



+l3^,Sf^,ap^T,{SfaP^CSf^C) }. (A7) 
For the common A, we have 

+7^75^7' o^pAC^Sf'^C'Ts^f T.cr''" + 7M75^r'7.^pATr(755fV^^C'5f^C) 
+7M75^r'^pAC^5r''^C'75^f 7.CT''" + 7M75^r'7.^pATr(755f a^^CSf "^C) 

+f5^,ST'a,^CSf^CSf^,aP^ + [i^.Sf ^,a,>:Ti{Sf a^^CSf^ C) }. (A8) 



3. Spin-3/2 Correlators 

Using the spin-3/2 interpolating field as given in Eq. ([T0|), the master formula is given 
for the nucleon: 

+2ST' lud^^^S'^' a^rCS::'^ Ccjpx) }. (A9) 

The master formulas for S and S are obtained in the same way as described above. 
For the octet A, we have 

r-p rp 

noaa' occ ohh', „aT i o Caa' occ ohh', „aT 

-Ib^ o^^Lb^ Lapxb.^ -iyO ^Ib^ a^r^b^ LUpxb^ ^^a 
+ASf ^^a'^^TiiSf a^rCS'f Capx) - Sf a,^CSf^Cap,St''^,a''^ 
+2Sf a^rCSf^CapxSf^^a'^- + Sf^^a'^^TiiSfa^^CSf^Cap,) 

naa' ^ ccd'^ „ nbb' ar n CO-ci' ^ Qcc'^ ^ cbb' aT 

-b^ a^r^Js ^(^pX^d ~ ^^u <^aT(-'Od ^<^pXJs If'^ 

+5r'7.^""Tr(5fV^.C5f }. (AlO) 
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For the singlet A, we have 

I T{ vt/.i^) ts,M }\^)- le'^^^e'^'^'^' a''^^ 

I Qaa' Qcc'^ ^ „ ebb', „ct qolo! „ Qcc'^ ^ Qbb', „<tt 

+Sf^,a''^TY{Sf'a,rCS':t'^Capx) }■ (All) 
For the common A, 

I T{ <=3/.(^) tU^) }\^)- le'^'^e^''' cj'^l.i 

-Sf- a^rCS"^ CapxS^^ ^^a"'' + S'^^ 7,,a''^Tr(5'f a^^CS^ Capx) 
-Sfa^rCSf^CapxSf-f^a^^ + Sf^^a^^TriSfa^rCSf^Capx) }. (A12) 
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